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Abstract 

The properties of nonlinear PDEs that generate filtered solutions are explored with particular 
attention given to the constraints on the residual term. The analysis is carried out for nonlinear 
PDEs with an emphasis on evolution problems recast on space-time-scale. We examine the role of 
approximation that allow for the generation of solutions on isolated scale slices. 

1 Introduction 

We are interested in applications, particularly of the evolution type, that are governed by nonlinear 
PDEs that generate solutions that exhibit fluctuations at all spatial/temporal scales. It is often the 
case that one is interested in obtaining a macroscopic solution that is defined as some kind of filtered 
representation of the dependent variables associated with some scale of resolution. Consequently, the 
macroscopic equations satisfied by the filtered variables will depend not only on space-time but also on 
a parameter associated with scale. We start by example and then generalise in the following sections. 

Let x = (x 1 , x n ) be the cartesian coordinates of 1™ (n = 1, 2, 3), representing the space domain, 
and denote the time coordinate by t S K+ (M + = the set of positive real numbers excluding 0). In order 
to define the filtered variables we introduce the scale parameter r\ G I on the scale interval / = (0,?7o) 
where < 770 << 1. We work on space-time-scale M = W 1 x K + x I C K m , m = n + 2, and introduce 
the linear operator 

£ =lr A (L1) 

where A = V • V and V is the spatial gradient operator. A filtered scalar variable 4> : M — » K can now 
be defined by 

C(j> = {) (x,t,r))€M. (1.2) 

Suppose that there exists i) (||$||l 2 (r™) < 00) such that lim^_ >0 + </> = $>(x,t). If we can also set 
lim^i^oo V(f> = then <f> is uniquely defined by 

<p(x,t,T])= I G(x-x,T])$>(x,t)dx 1 ...dx n (1.3) 



where G(x, rj) = {\itr\)~ n l' 2 exp[— |a;| 2 /(4?y)] with the property J Rn G(x, if) dx 1 . . . dx n = 1. Hence 4>{x, t, rf) 
becomes the Gaussian filter of &(x, t). In this way we can identify r/ — /3S 2 , where S is a parameter 
associated with the characteristic spatial resolution and j3 is a constant associated with the filter width. 

A major sticking point with this approach is the need to assume the existence of a regular limiting 
solution, as r\ — ► + , associated with the space-time pointwise representation of the dependent variables. 
Here we shall depart from the approach made in [3] and [3] by attempting to work entirely in M and 
avoid any reference to a limiting fully resolved solution on space-time. 

We denote by 1Z(M, 1^), N > 1, the set of smooth vector functions on M, i.e. smooth mappings 
M -> R N . For the case N = 1 members of 1Z(M, M. N ) are smooth scalar functions. 

Consider now an ideal fluid whose variables are defined on M as follows. Let the fluid velocity 
v e K(M,R n ) be a filtered field, defined by 

£v = (x,t,r])eM, (1.4) 

and satisfy the continuity condition 

V-u = (x,t,r])eM. (1.5) 
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We now construct the macroscopic equations from (|1.4p and (|1.5p . 

On each scale slice M| 7)=const , in a time interval Icla fluid element follows a path under the flow 
characterized by a mapping 7 : 1 — > M\ v=const (integral curve of v) so that 

dx o j(t) 

vo ^ t ) = — d f J -- ( L6 ) 
The acceleration, a, of the same fluid element can be defined along the flow path j(t) by 

ao 7(i ) = ^W=^). (1.7) 
Alternatively we may write a € 1Z(M, 1") as 

a = — +v ■ Vv. (1.8) 
at 

We introduce the symmetric second order tensor a whose components are given by 

a a b = f2^1^ l<aih<n . (L9) 

^ dx c dx c - ' - v ; 

c— 1 

We also define 

s = -2V-cr. (1.10) 
The properties assigned to v above dictate that the acceleration must satisfy the following constraint. 

Proposition 1.1 If v e 1Z(M, R n ) satisfies fTfl ) and fT5j) then 

Ca = s (x,t,r))€M, (1.11) 

where s is given by 11.10\) . 
Proof. We have 

Ca = + v ■ Vv) = —(Cv) + Civ ■ Vv) . (1.12) 

at at 

By a straightforward calculation 

£(v ■ Vv) = {Cv) ■ Vv + v ■ V(£«) + s (1.13) 

where we make use of the fact that v is divergent free to express the source term s in the divergent form 
(|1.10[) . The required result follows from (|1.12[) and (j!.13|) and the fact that v is a filtered field . | 

We can always introduce a filtered scalar function p G 1Z(M, R), 

£p = (a;,t,?j)eM, (1.14) 

and write 

a = -Vp + r, (1.15) 

for some r £ 1Z(M, R n ). Noting that £a = —CVp + Cr = —VCp + Cr = Cr, we have constructed from 
(|1.4[> and (jl.5|l the system of macroscopic equations: 

V-v = {x,t,n)eM, (1.16) 

-7^ + v • Vv + Vp = r {x,t,n)eM 7 (1.17) 
£r = s (x,t,r))€M, (1.18) 

where s is given by (ll.lOp . 

We can think of the residual r as incorporating the effects of the large scale viscous and stress/strains 
of the fluid that are manifestations of the filtering process. In the case that we can assume the existence 
of a limiting pointwise representation of the filtered variables we can impose, under some appropriate 
norm || • ||, the condition that lim I) _ ! .Q+ \\r — isAv\\ = 0, for some constant v. In this way we recover, at 
least in the generalised sense, the Navier-Stokes equation from the macroscopic equations (|1.16p - (|1.17p 
as a limiting formulation. Since the notion of viscosity is itself meaningful only in some locally averaged 
sense, under the continuum fluid model the inviscid limit, lim r/ ^ + IMI = 0, is more natural. 
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2 Macroscopic Equations 



We now attempt to generalise avoiding application specific issues. The system of nonlinear PDEs con- 
sidered will be of first order but it should be apparent that the analysis can incorporate higher order 
systems. As in the previous section, let x — (x , x n ) be the cartesian coordinates of M n (n — 1, 2, 3), 
representing the space domain, and denote the time coordinate by t G R+ (M+ = the set of positive real 
numbers excluding 0). Let M = M. n x K + x I c R"\ m = n + 2, where I — (0, 770) is the scale interval, 
77 G I is the scale parameter and < 770 << 1. We shall write £ = (x,t,ri) — . . . , £ m ), where we have 
set £ a = x a (1 < a < n), = t and £ m = £"+ 2 = r\ is the scale parameter. 

For any D C R^ 1 , > 0, we denote by TZ(D,R N2 ), N 2 > 0, the set of smooth (C°°) maps it : D ->• 
I™ 2 . For it = = (it 1 , ...,1^) G 7?.(M, R N ) we shall denote u z = d % u = «), u tJ = d l d j u = (ug), 
. . . , with the shorthand notation di = d /<9£\ Throughout, unless otherwise stated, we use the range 
< i,j,k < m for lower case Latin indices and the range 1 < a, /3,7 < TV for the lower case Greek 
indices. For easier identification of certain important operators we use the indices t instead of n + 1 and 
r\ instead of m = n + 2, e.g. we write dt instead of d n +i and d v instead of d n +2 — 9 m . 

We shall express our macroscopic equations in terms of a core function of u and Ui that is form 
invariant with respect to scale and a residual term. Let F e M. N such that F(u, ui) is a smooth vector 
function with respect to its indicated arguments. Given any u £ 1Z(M, l w ) there exists an r S 1Z(M, R N ) 
such that 

F(u,u z )=r ££M. (2.1) 

Alternatively, we may regard (|2.1[) as a system of N first order PDEs for the unknown dependent vector 
function u(£) if r(£) is known. We shall henceforth refer to r as the residual and F(u,Ui) as the core 
function. The core function may be regarded as retaining the structure of the N PDEs that are associated 
with the fully resolved system as shall be outlined below. The construction of the core function in a 
wider context is application specific and follows along the lines discussed in the previous section. We 
shall return to this application in a later section. Our first task is to find constraints on r(£) such that 
the solutions of (|2.1[) have well defined filter properties. 
We introduce the vector field operators Vi defined by 

Vi=di + ufd a + ufjdi + u% k 8i k + ... (2.2) 

with the notation 

d ■ d • d • ■ ■ d 

dl = W da ^d^' 9 " = a<' 5 " = 9<' "' (2 ' 3) 

Here and throughout we use the summation convention of repeated upper and lower indices over the 
assigned ranges for the given indice type given above. For u G 1Z(M,M. N ) the operator Vi acting on 
any F G M. N using the more general representation F(£,u,Ui,Uij, . . .) generalises the partial derivative 
operator di acting on members of TZ(M, Mr). {V%F is sometimes referred to as the total derivative of F 
with respect to The action of Vi on any w G TZ(M, M N ), w — w(£), gives the partial derivative of w 
with respect £*, i.e. V,w = diW. The vector field operators Vi may be truncated up to terms associated 
with the largest order of the partial derivatives of u present in the expressions on which Vi is operating. 
In truncated form the vector fields Vi have special geometric significance in the general theory of PDEs 

mm- 

We introduce the spatial Laplacian operator 

A = V-V, (2.4) 

where V = (<9i,-- - , d n ) is the spatial gradient operator. The associated operator acting on vector 
functions F G Mr using the general representation F(£, u, u%, Uij, . . . ) is defined by 

n 

L = J2 VbVb ■ ( 2 - 5 ) 

6=1 

The set of filter maps is defined as follows: 

Definition. P(M,R N ) is the set of vector functions u G TZ{M,R N ) such that 

(dn-A)u = £eM. (2.6) 
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We can assume that P(M,~R N ) is sufficiently populated that it is of interest to the analysis that 
follows. Suppose for the moment that we can assume that there exists u £ 7\L(R" x R+jR^) satisfying 
the first order PDEs 

F(u, d l u) = (x, t) £ R™ x R+ (2.7) 

for some F £ W N such that F(u, ui) is smooth with respect to its indicated arguments. If for some u £ 
V(M, K) we have u\ v=0 + = u and limui^oo Vu = then u(x, t, rf) = L„ G(x — x, r))u(x, t) dx 1 . . . dx n , 
where G{x,rf) — (4nri)~ n / 2 exp[— |x| 2 /(4?7)]. Hence it(£) — u{x,t,rf) is the Gaussian spatial filter of 
u(x, t). 

The PDEs (|2.7p represent the fully resolved system. It should be mentioned at this point that in 
many applications of interest F will not explicitly depend on terms u v (see application below). Under 
our assigned range for the lower case Latin indices the representation F(u,Ui) will indicate a possible 
explicit dependence on u v . The more general representation F(u,Ui) will not effect the calculations 
below and is used to avoid introducing more notation for the indices and their ranges. 

Since for u £ V(M, M. N ) we cannot expect that in general F(u,Ui) will vanish identically on M we 
must introduce a residual term r € 1Z(M, R w ) as in (|2.ip . The system (|2.7p is recovered from (|2.1[) by 
requiring that lim, ) _ s .g+ r = 0. 

In many applications one cannot gaurantee the existence of regular solutions u € lZ(M. n x R+jR^) 
for the fully resolved system and it is at this point that we depart from the assumptions made in [3] , [4] . 
Our focus is on nonlinear problems that generate solutions that exhibit fluctuations at all space and time 
scales. We propose that such systems can be defined in a meaningful way enirely within space-time-scale 
avoiding any reference to a fully resolved system. 

We shall make extensive use of the vector field operator 

W = d n + {Au a )d a + (Auf)fl* + [AufyfrJ + ... (2.8) 

We are now in a position to establish the constraints on the residual for filter maps. 

Proposition 2.1 If u e V(M,R N ) and r G H(M,R N ) such that F(u,Ui) - r = 0, £ e M, where 
F G M. N and F(u,Ui) is smooth with respect to its indicated arguments, then (d v — A)r = s, £ G M, 
where s = (W — L)F(u, Uj). 

Proof. For any u £ 1Z(M, M. N ) there will always exist an r £ 1Z(M, M. N ) such that F(u, Ui) 
identically on M. We have 

(V v -L)(F(u, Ui )-r)=0. 
Manipulating the left hand side we obtain 

(V v - L){F{u,u t ) -r) = (V v - L)F(u, Ui ) - {8 V - A)r 

= (V n - W)F(u, Ui ) + (W- L)F(u, Ui ) - {8 V - A)r (2.10) 
= % - W)F(u, Ul ) + s- {d v ~ A)r , 

where we take r £ 1Z(M,M. N ) to mean that it has the representation r(£) so that Vtr — dir. It is 
straightforwad to see that if u £ V(M,M. N ) then V v — W and hence the first term in this last expression 
of (|2~T0|) vanishes. Thus 

(V v - L){F{u, Ui) -r) = -(d v - A)r + s . (2.11) 
The desired result follows from (j2~9| and (j2~TTj) . | 

We can now write the full system of macroscopic equations for the filtered variables u £ V(M,M. N ) 

as 

F{u,u t )=r £eM, (2.12) 
(d v -A)r=s $eM, (2.13) 

where 

s = (W - L)F(u, u.^ . (2.14) 

Given any core function F £ R N such that F(u,Ui) is smooth with respect to its arguments, for 
any u £ V(M, M N ) there will always exist an r £ K(M,R N ) in (|2~T2l) satisfying ([2T3]) - ([2T4ll . From the 



— r vanishes 
(2.9) 
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practical viewpoint, we must consider as a separate issue the task of developing algorithms that extract 
from V(M, Mr) those members that can be regarded as meaningful, in some formal sense, to the specific 
application under consideration. 

We now establish an important link between the filter error, if), and the residual equation error, e, 
as defined in the following. 

Proposition 2.2 Let u, r £ H(M, M N ) such that F(u, m) - r = 0, £ £ M, where Fel" and F{u, m) 
is smooth with respect to its indicated arguments. Let ip £ 1Z(M, R ) be defined by 

i/>(£) = (d n - A)«(0 . (2.15) 
It follows that ip satisfies the linear system of first order PDEs 

Ctfditf + =e a £ € M, (2.16) 

where the coefficients are given by CT = d l pF a (u 1 Uj), — dpF a (u,Ui), and 

e=(d v -A)r-s, (2.17) 

where s — (W — L)F(u, Ui). 

Proof. We follow a similar argument to that of the proof of the previous proposition but note that now 
r need not satisfy (|2.13p - (|2.14p and hence u is not necessarily a member of V(M, M. N ). Since F(u, Ui) — r 
vanishes identically on M we have 



(V v -L)(F(u, Ui )-r)=Q. (2.18) 



The left hand side becomes 



{V v - L){F(u, Ul ) ~r) = (V v - L)F{u, Ui ) - (d v - A)r 

= (y v - W)F(u, Ui ) + (W- L)F(u, m) - (d v - A)r 
= (V v - W)F{u, Ui ) + s- (d v - A)r 
= (V v - W)F(u, Ui)-e, 



(2.19) 



where we take r £ 1Z(M, K> N ) to mean that it has the representation r(£) so that Vir — dir. The first 
term in the last expression will not in general vanish because u is not necessarily a member of V(M, Mr). 
Using the definitions of the vector field operators V v and W, we have the explicit representation 

V v - W = ^ a d a + (d^ a )0 l a + (did^dH + ... (2.20) 
The identity ([2~16| then follows from pT8]l - ([230]l . | 

If the residual r satisfies ()2.13|) . where the source term is given by (|2.14j) . then e vanishes identically 
on M and ()2.16j) admits the trivial solution from which it follows that u £ V(M, M ). 

The system of PDEs (I2.12p ~ d2.13p are the exact equations satisfied by the filtered variables u £ 
V(M, M ) on space-time-scale M. They are of little use in applications since the presence of the d n 
term in the residual equation (|2.13|) means that the solution of the system can only be obtained by scale 
integration from some arbitrary scale slice of M on which u is prescribed. To overcome this problem we 
resort to approximation by replacing the residual equation (|2 . 1 3|) with a constraint on the residual that 
allows us to generate regular solutions u independently on scale slices M| j;=cons t such that u is a good 
approximation of some member of V(M, M ) on that scale slice. 



3 Approximation 

An example of a residual approximation that has been found to be useful in applications is one in which 
the exact residual equation (|2.13l) is replaced by the constraint (see [3],[1]) 

Ar-- + s = 0, (3.1) 
V 
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where the source term s is given by (|2.14l) . The approximate macroscopic system (|2.12l) and (|3.1[) can 
be solved independently on any scale slice M|^ =const by updating the residual via (|3 . 1 1) at each timestep 
as soon as u is available. If we can assume that lim 

dr r 
drj r\ 



r)->0+ 



V 

z M 



it follows that 



d 2 r 



drf 



(3.2) 



Motivated by a discussion of residual approximations in general, we shall avoid any focus on specific 
formulations such as fl3.1[) . 

In practical application space-time solutions are generated from the system of PDEs (|2.12l) on some 
scale slice M| t?=const while using a residual that is not exact, i.e. one in which e does not vanish identically 
on that scale slice. If we enforce (|2 . 1 2[) with a residual that is not exact then it follows from Proposition 
12.11 that u cannot be a member of V(M, R ). 

Our objective here is to investigate how well members of 1Z(M, M. N ) serve as approximations to mem- 
bers of V(M, M. ). More precisely, suppose that we employ a residual that is based on an approximation 
of the residual equation (|2.13l) for which we know that a regular solution u G 1Z(M, ~R N ) for (|2 . 1 2[) 
exists. We are interested in how u will deviate from some member of V(M, M. N ) as the scale parameter 
is increased if the two solutions agree on some arbitrary scale slice. 



Proposition 3.1 For some e G (0,%), let M e = M\ ve{e , Vo) . Let u G K(M e ,R N ) and u G V(M, R") 
such that lim ?7 _ >e + (u — u)=0. Suppose further that tp G 1Z(M £ ,M. N ), defined by = (d v — A)u(^), is 
bounded on M s and that limij-i^oo Vu, Vu = 0. It follows that 

(u — u)[x,t,rj) = / / G(x — x, rj — fj)tp(x, t, fj) dx 1 ... dx n di) (x,t,rj) € M s , (3-3) 
Je JR n 

where G(x,T)) = (47rn)~ n / 2 exp[—\x\ 2 / (4n)]. 



Proof. Setting cj = u — u we have 



We also have 



lim Vw = . 



(3.4) 
(3.5) 



The solution of (I3.4p - ()3.5|) is unique and is given by p.3[) . | 

Given that J R „ G(x, rf) dx 1 . . . dx n = 1 we obtain from p.3p the bound, for ^ G M e 

\u — u\ < rj sup \ip\ . 



(3.6) 



4 Application 

We return to the example of an ideal fluid. We set u — (v,p), where v — (v , . . . , v n ) is the filtered 
fluid velocity vector and u n+1 = p is the filtered fluid pressure. In this case N = n + 1 and we set 
F = (F ( „),F (p) ), where F (v) = (F 1 , . . . , F n ) and F (p) = F n+1 . Motivated by Section 1 we define the 
core function by 

F {v) = v t + v-Vv + Vp, F {v) =V-v. (4.1) 

Note that the core function is not explicitly dependent on £ and u^. 

Since we insist that admissible solutions have well defined filter properties, i.e. u G P(M,M. N ), 
Proposition 12 . II demands that the residual satisfy (|2.13[) . Here the source term s is obtained by applying 
the identity (12. 14)) to the explicit representation of the core function given in (14. 1[) : 

n n 

sa = - 2 EE^<' (l<o<n), S " +1 =0, (4.2) 

fc=l c=l 



G 



using the notation u, = («i,J>i), Uij = (vij,Pij). We set r — (rr v ) , J7 p \ ) , where rv„) = (r , • • • , r n ) and 
r (p) = r " +1 : an d similarly s = (s(„), S( p )), where S(„) = (s , • • • , s n ) and S( p ) = s n+ . If we can assume 
that rVpj — >• as i] — > + then, in light of the residual equation (|2.13[) and the fact that the component 
of the source term s/ p ) = s n+1 vanishes identically on M, ru,\ also vanishes identically on M. It follows 
that the continuity equation V ■ v = is form invariant with repsect to scale. We can now write 

s (w ) =-2V-(r, S( P) =0, (4.3) 

where a is a symmetric second order tensor whose components are given by (jl.9l) . Thus (|4.3p is in 
agreement with the results presented in the first section. We note that the source term in the residual 
equation also does not contain any explicit dependence on terms involving d n and dt- 
For u £ 'P(M, M N ) the macroscopic equations (|2.12p then take the explicit form 

dv 

— +v-Vv + Vp = r [v) , V • v = . (4.4) 

In practice we wish to solve (j4.4[) (subject to suitable initial conditions) independently on scale slices 
M\ n=const . Although the source term (|4.3[) is suitably defined, the residual equation (|2 . 1 3|) is still not 
suitable for this purpose due to the presence of the d v term on the left hand side. For this reason we 
resort to an approximation of the residual based on ([2.13|) - ()2.14|) with the aim of obtaining a member 
of H(M,M. N ) that approximates some member of V(M,R N ). 

In a similarly fashion as above, let us set e = (er v ),e(p)) and tp = (ip(v) ? V'(p))- Because rv p ) and S( p ) 
vanish indentically on M it follows that ei p \ also vanishes identically on M. Using the prescription of 
the core function (|4.1j) to obtain the coefficients, the system of PDEs (|2.16j) reduce to 

dtip{ v ) + v ■ V^ ( „) + t/»(«) • Vw + V-0( P ) = &{ V ) , V • ip( v ) = . (4.5) 

for some v £ 1Z(M, K") that is a solution of (|4.4j) . Clearly the long time behaviour of (|4.5|) is of interest 
and will depend crucially on the specific residual approximation used. If one can supply reliable estimates 
of the residual equation error e on the scale slice on which the space-time solution is being sought then 
there is the opportunity of resorting to computation. Since the left hand sides of (|4.5p contain no partial 
derivatives with respect to scale we may, from the computational view, obtain numerical estimates of ip 
by coupling (|4.5|) with the system (|4.4p using the same algorithms. In this way one may obtain estimates 
of the error bound (|3.6p on the scale slice of interest as the computations proceed. 

Before concluding it is important to note that for any u £ V(M,'M. N ) satisfying the continuity 
constraint there will always exist an ri v -s £ V(M,M. N ) satisfying (12 . 13[) with source f)4.3[) such that (|4.4p 
will hold. This does not entirely resolve the existence problem for applications since many solutions 
will have very little interest in this respect. As already mentioned, the problem of extracting from this 
set those solutions that are of interest to the application under study is a separate issue and must be 
examined in the context of the prescribed initial conditions and, to some extent, the solution method. 

5 Concluding Remarks 

In avoiding reference to a fully resolved solution we do not insist that regular fully resolved solutions do 
not exist, only that we do not insist that our space-time-scale solutions uniformly converge to one in the 
limit as rj — > + . In doing so we have the expectation that meaningful solutions to applications may be 
obtained entirely within space-time-scale. 

We have already seen that for any member of V(M, R w ) there will exist an r £ 1Z(M, R N ) satisfying 
(|2.12p and (|2 . 13[) where the source term s is given by (12 . 14[) . But we have yet to formally establish that 
V(M,M. N ) is sufficiently populated that it contains members that have any useful meaning for specific 
applications. Allied to this is the construction of algorithms that allow us to extract from V(M, M. N ) 
such solutions. 

Some progress has been made in this direction through numerical solution methods. For nonlinear 
systems the introduction of discretisation dictates that one is seeking an approximation of the filtered 
representation of the dependent variables associated with the scale of resolution rather than an approx- 
imation of the fully resolved dependent variables. In this context there arises the associated problem of 
how one defines the core function. We can no longer simply associate the core function with the PDEs 
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that are to be assumed to govern the fully resolved system. In Section 1 an attempt was made to redcrivc 
the core function of an ideal fluid application without reliance on the well known governing equations 
assumed for the fully resolved system. At this stage, the problem of defining the core function is very 
much application specific and needs to be investigated further in a more general context. 
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